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ABSTRACT 

Although the gravitational wave kick velocity in the orbital plane of coalescing 
black holes has been understood for some time, apparently conflicting formulae 
have been proposed for the dominant out-of-plane kick, each a good fit to different 
data sets. This is important to resolve because it is only the out-of-plane kicks 
that can reach more than 500 km s _1 and can thus eject merged remnants from 
galaxies. Using a different ansatz for the out-of-plane kick, we show that we 
can fit almost all existing data to better than 5%. This is good enough for any 
astrophysical calculation, and shows that the previous apparent conflict was only 
because the two data sets explored different aspects of the kick parameter space. 

Subject headings: black hole physics — galaxies: nuclei — gravitational waves 



1. Introduction 



When two black holes spiral together and merge, the gravitational radiation they emit 



to the original center of mass. The speed can reach over 3,000 km s 
which would eject the remnant from any galaxy in the universe (see Figure 2 of lMerritt et al. 



Dain et al. 


2008 




s2of 


Merritt et al. 



1 CRESST and Gravitational Astrophysics Laboratory, NASA/GSFC, Greenbelt, MD 20771 

2 Department of Physics, University of Maryland, Baltimore County, 1000 Hilltop Circle, Baltimore, MD 
21250 

3 Department of Astronomy, University of Maryland, College Park, MD 20742 
4 Joint Space Science Institute, University of Maryland, College Park, MD 20742 
5 Gravitational Astrophysics Laboratory, NASA/GSFC, Greenbelt, MD 20771 
department of Physics, University of Maryland, College Park, MD 20742 



-2- 



( 12004 ) ). As discussed in iMerritt et al.l (12004 ). the kick magnitude and distribution are im- 
portant for discussions of hierarchical merging, supermassive black hole formation, galactic 
nuclear dynamics, and the degree to which black holes influence galaxy formation. The 
community has converged on the formula f or the kick speed in the original orbital plane 
( IBaker et al.l 120071 ; ICampanelli et al.l l2007bl ; I Gonzalez et al.l 120071 ) , but apparently conflict- 
ing dependences for the out-o f-plane kick (which d o minat es the total kick for most configu- 
rations) have been proposed. lLousto fc Zlochowerl (120091 ) suggested that for a binary with 



component masses m x and m 2 > m 1; the kick sc ales as r) 2 , w h ere rj = m 1 m 2 /(m 1 + m 2 



is the symmetric mass ratio; however the data in IBaker et al.l (120081 ) were fit much better 
with an if dependence. The conflict is only apparent, however, because there were no runs 
in common between the two data sets. This suggests that an analysis might be performed 
with a new ansatz that can fit all of the existing data. 

Here we perform such a fit, and demonstrate that there is a single formula for the out- 
of-plane kick that fits almost all exist ing data to better th an 5% accuracy. The new ansatz 
is similar to one recently suggested by lLousto et al.l ( 120091 ). for which, however, a fit was not 
attempted. Its form is based straightforwardly on the post-Newtonian (PN) approximation 
and includes both the aforementioned if and if terms, as well as a slightly more complicated 
spin-angle dependence. In § 2 we describe our new runs and we describe the new ansatz in 
§3. In § 4 we list all 95 runs we have fit, from different numerical relativity groups, and our 
fitting procedure and best-fit parameters. In § 5 we discuss the implications of our results, 
and indicate the fraction of kicks above 500 km s _1 and 1,000 km s _1 for representative spin 
and mass ratio distributions, comparing it with previous results. The goodness of our fit 
to the entire usable data set of out-of-plane kicks suggests that the full three-dimensional 
kick is now modelled well enough that it will not limit the accuracy of any astrophysical 
calculation. 



2. Numerical simulations 



We have performed new simulations representing 22 distinct physical cases. Defining 
q = mi/ m2 and ojj = Si /in 2 where Si is the spin angular momentum of the ith black hole, 
we used mass ratios of q = 0.674 or q = 0.515, and spins initially within the orbital plane 
of a% — (%i — 0.367 and a 2 = = 0.177 or a 2 = ctf" = 0.236, where the _L superscript 
indicates orthogonality to the orbital axis. For each mass ratio we used one of eleven different 
spin orientations (given in Table 1), in order to probe the spin-angle dependence of the final 
recoil. Initial parameters were informed by a quasicirc ular PN approximat ion and initial data 
constructed using the spectral solver TwoPunctures ( lAnsorg et al.l 12004 ). Evolutions were 
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performed with t he Einstein-solver Hahndol ( llmbiriba et al.l I2004J : Ivan Meter et al.l 12006 
Baker et al.l 120071 ). with all finite differencing and interpolation at least fifth-order-accurate 
in computational grid spacing. Note for the purpose of characterizing the simulations, it 
is convenient to use units in which G — c — 1 and specify distance and time in terms of 
M = mi + m2- The radiation field represented by the Weyl scalar \l/4 was interpolated to 
an extraction sphere at coordinate radi us r = 50M and integr ated to obtain the radiated 
momentum using the standard formula (ISchnittman et al.l 120081 ): 

2 



P; 



dt- 



16vr 



dQ 



X ■■>) 



dtm, 



For each of the 22 cases we ran two resolutions, with fine-grid spacings of hf — 3M/160 
and hf = M/64. To give some indication of our numerical error, the total recoil from the 
two resolutions for each of the q = 0.674 cases agreed to within ,< 6%, and for each of the 
q = 0.515 cases to within 1% (the latter using a more optimal grid structure). 

In these simulations, variation in the magnitude of the final recoil within the x-y plane 
suggested non-negligible precession of the orbital plane. Indeed the coordinate trajectories 
of the black holes showed precession of up to ~ 10°. To calculate the component of the recoil 
velocity parallel to the "final" orbital axis, Vu , we tried two different methods. In one method 
we simply took the dot product of the final, numerically computed recoil velocity V with the 
normalized orbital angular momentum L|L| _1 as calculated from the coordinate trajectories 
of the black holes, just when the common apparent horizon was found: V\\ = V ■ L|L| . In 
our second method we assumed that each black hole spin, which is initially orthogonal to the 
orbital angular momentum, remained approximately orthogonal throughout the simulation, 
i.e. 



L ■ Si ~ L ■ 



0. 



(2) 



This assump tion is consistent with PN calculations, to linear order in spin (jRacinel 12008 



Kidderl Il995l ). and is also supported numerically by the fact that the merger times in our 
simulations were independent of the initial spin, to within At < 1M. In this case the in-plane 
recoil should depend onl y on the mass r a tio (a nd not the spin), and we assume it is given 
by the formula found by I Gonzalez et al.l (120071 ) , with coefficient values given by a previous 
fit (IBaker et al.l 120081 ). This implies 



V 2 - {V_ 



pred\2 



IVJ 



(3) 



where we have further assumed that the sign of V\\ should be the same as that of V z , given 
the modest amount of precession. 

These two definitions for Vu were found to differ by a relative error of <, 5% for all points 
except one for which Vu V±. Relative to the maximum Vu per mass ratio, they were found 
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to differ only by < 2%. Note this lends further support to our assumption that the spins 
are orthogonal to the orbital angular momentum to good approximation throughout these 
simulations (Eq. (T5])). Even V z was found to differ from the above definitions for V\\ by only 
< 2%, relative to the maximum. We will use Eq. (jHJ) to give our canonical V\\, for the purpose 
of analytic fits to the data. 

An additional assumption we will make about our data, important for fitting purposes, 
is that the amount of precession undergone by the spins in the orbital plane, i.e. the dif- 
ference in spin-angles, between the initial data and the merger, is independent of the initial 
spin orientation. This assumption, valid to linear order in spin according to the PN approx- 
imation (Eq. (2.4) of (lKidderlll995l )). was previously fou nd to be the case t hrough explicit 
computation of the spins in the simulations presented in iBaker et al.l ( 120081 ). We have not 
explicitly computed the spins in the new simulations presented here but the independence 
of the merger time with respect to initial spin orientation is consistent with the assumption 
of similar independence of spin precession. 

For the purpose of constructing an accurate phenomenological model, we added to this 
data set previously published data representing out-of-plane kicks. Our criteria for selecting 
relevant data are that the component of the final recoil parallel to the orbital plane is given 
and at least three different in-plane spin orientations are used. We arrived at a t otal o f 
95 data points: the 22 new ones presented here plus 73 drawn f rom IBaker et al.l ( 120081 ); 
Campanelli et al.l d2007bh: bain et all J2OO8I ): lLousto fc Zlochowerl J2OO9J). Note that for the 



cases m 



Lousto &: Zlochowerl ( 120091 ) we use the second version of their values listed in their 



Table III, resulting from what they consider to be their best calculation of the final orbital 
plane. 



Ansatz 



An ansat z for the total reco i l that was found to be very consistent with num erical results 
has the form ( IBaker et al.l 120071 ; ICampanelli et al.ll2007bt iGonzalez et al.l 120071 ) 



V 



Vxm e i + Vj_ s (cos£ei + sin£e 2 ) + V\\ e 3 , 



477(1 



H- 



1 + q) 



Q. 



(4) 
(5) 

(6) 



where ei and e3 are unit vectors in the directions of separation and the orbital axis just 
before merger, respectively, e 2 = ei x e3, a\ and a| represent the components of spins 
parallel to the orbital axis, £, A, B, and H are constant fitting parameters, and V11 is to be 
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discussed. 



Every component of Eq. fl5j) has been modeled to some degree after PN expressions. 
The use of such PN-based formulae has proven very successful. For example, Eq. (J5J) for 
the in-plane, mass-ra tio-determined r ecoil can be obtained from the corresponding PN ex- 
pression (Eq. (23) of iBlanchet et al.l (120051 )) simply by taking the leading-order terms and 
replacing instances of the PN expansion parameter (in this case frequency) with constant 
fitting parameters. Using this. iGonzalez et al.l (120071) obtained very good agreement with a 
large set of numerical results. Le Tiec et al. ( 201oh . who calculated the recoil using a com- 
bination of the PN method with a perturbative approximation of the ringdown, also found 
that Eq. (jSJ) gave a good phenomenological fit to their analytic results. 

Why such a prescription for generating an ansatz from the PN approximation should 
apply so well to merger dynamics is not perfectly understood. The effective replacement 
of powers of the frequency with constants may be defensable because the majority of the 
recoil is generated within a narrow time- window near merger, perhaps within a narrow range 
of frequencies. This is particularly evident for the out-of-plane recoil speed, which rapidly 
and monotonically increases up to a constant value around merger. However, values of 
the constant coefficients that appear in the PN expansion cannot be expected to remain 
unchanged because, as merger is approached, high-order PN terms with the same functional 
dependence on mass and spin as leading order terms can become comparable in magnitude. 
So, it is just for this functional dependence that we look to the PN approximation for 
guidance. 

Followin g the prescrip ti on im plied above, from the PN expression for out-of-plane re- 
coil given by iRacine et al.l (120091 ). Eqs. (4.40-4.42) (neglecting spin-spin interaction), the 
following can be straightforwardly obtained: 

K 2 r] 2 + K 3 r] 3 



Vn 



+ 



g+1 
K s (q ~ l)if 



[qoci cos( 



- $i) - or 2 



cos 



-*2>] 



2 a^ cos( 



— $i) + ak cos 



(7) 



where K 2 , K 3 , and K$ are constants, af- represents the magnitude of the projection of the 
ith black hole's spin (divided by the square of the black hole's mass) into the orbital plane, 
(pi represents the angle of the same projection, as measured at some point before merger, 
with respect to a reference angle representing the direction of separation of the black holes, 
and represents the amount by which this angle precesses before merger, which depends 
on the mass ratio and the initial separation. We have ignored terms quadratic in the spin 
be cause we assume the y are subleading. Note that this ansatz is equivalent to one suggested 
by lLousto et al.l (120091 ) . provided the angular parameters are suitably interpreted. In terms 
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of the notation of iRacine et al.l (120091 ) . 



(5 + 1) 1 [q®i cos 
(g + 1)" 2 [q 



a 1 cos^i 



$1) 



a 2 cos( 

Ct COsf 



-$2)] 
-$2)] 



-M- X A ■ n 
M" 2 S • n 



-M- 1 A ± cos(6),(8) 
M- 2 ^ ± cos(^), (9) 



where n is a unit separation vector, A 
l^K - Q<4\, S L = M 2 (q + l)~ 2 \ai 4 



2 



S 2 /m 2 - Si/m 1} S = Si + S 2 , A x = Af(g + 
is the angle between A and n, and \I/ is 



the angle between S and n, all measured, for our purposes, at some point arbitrarily close to 
merger. It may be interesting to note that the expression multiplying K$ takes into account 
effects of orbital precession (because it is non-vanishing if and only if the orbit precesses) 
and therefore represents physical phenomena neglected by previous fits. 



4. Fitting procedure and results 

The fitting of the out-of-plane recoil is in principle complicated because in addition 
to the overall factors K 2 , K 3 , and K$ (which are the same for any mass ratio or spins), 
any particular set of runs with the same initial separation and mass ratio (which we will 
term a "block") has idiosyncratic values of $1 and $2 that, although not fundamentally 
interesting, need to be fit to the data. Therefore, in the 17 blocks of data we fit, there are 
formally 3 + 2 x 17 = 37 fit parameters. In the eight data blocks for which = 0, $1 never 
enters, and in the two for which a± = a 2 , $1 = $2- Therefore the actual number of fitting 
parameters is 27, but this is still large enough that a multiparameter fit would be challenging. 
Fortunately, each ($1, $ 2 ) pah only affects a single data block. We can therefore speed up 
the fitting, and incidentally concentrate on only the interesting parameters, if we (1) pick 
some values of K 2 , K 3 , and K$ that apply to all data blocks, then (2) for each data block, 
find the values of $ 2 and possibly $1 that optimize the fit, repeating this using new values 
of $1 and $ 2 for each block. This gives an overall fit for the assumed values of K 2 , K3, and 
Kg, having optimized over the uninteresting <£> parameters. 

The fit itself needs to be performed assuming uncertainties on each of the numerical 
measurements of the kick. Each such calculation is computationally expensive and systematic 
errors are usually difficult to quantify, hence we do not have enough information to do a true 
fit. As a substitute, we assume that for each block of data, the uncertainty a in each kick is 
either equal to a fraction (fixed for all blocks) of the maximum magnitude kick in the block, 
or to a fixed fraction of the individual kick itself. The former may be justified because some 
sources of error will be independent of the phase of the angles when the holes merge, but 
we note that the fit performed with the latter assumption (that the uncertainty equals a 
fractional error of each kick) yields very similar values for the fitting factors. As we do not 
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know what the actual fractional error is, in either case we adjust it so that for our best fit 
we get a reduced x 2 °f roughly unity (given our 95 data points and 27 fitting variables, this 
means we need a total \ 2 °f about 68). We then evaluate every (K 2 , K 3 , Kg) triplet using 
X 2 = XXpred — kick) 2 /cr 2 . 

Minimizing \ 2 as calculated with respect to the maximum kick per block, we find that 
a 2 = 0.0005V^ 2 max (block) gives x 2 /dof = 1.0. We note that this value for a is comparable to 
the numerical error as measured by the difference in kicks computed at different resolutions, 
when available (e.g. for the new simulatio ns presented here, or the q = 0.25 case presented 



in Table IV of lLousto fe Zlochowerl ( 120091 )). It is also worth noting that for the data from 
the new simulations, the uncertainty due to orbital precession, discussed in § [2], is less 
than 0.02Vji max < a, and therefore cannot significantly affect the fit. Our best fit is K 2 = 
30,540 km s~\ K 3 = 115,800 km s _1 , and K s = 17,560 km s -1 , with la ranges 28,900 - 
32,550 km s~\ 107,300 - 121,900 km s _1 , and 15,900 - 19,000 km s" 1 , respectively. We 
emphasize that all three coefficients are indispensable in obtaining a good fit; e.g. using the 
same definition of \ 2 as above, if K 3 = then the best fit gives x 2 /dof = 2.0, or if Ks = 
then the best fit gives x 2 /dof = 1.6. 

Minimizing x 2 as calculated with respect to each individual kick, we obtain similar 
results. In this case, a 2 = 0.0016 V|p gives x 2 /dof = 1.0. Our best fit becomes K 2 = 
32,092 km s _1 , K 3 = 108,897 km s -1 , and Ks = 15,375 km s _1 , in agreement with the 
above fit to within ~ 5% for K 2 and K 3 and < 12% for Ks- 

In Table H] we compare the predicted out-of-plane kicks with the measured ones for the 
entire data set. For fit #1 minimizing the error with respect to the maximum kick per 
block, of the 95 points, only 4 agree to worse than 10%, and all of those occur for kicks 
with magnitudes much less than the maximum in their data block. Therefore these could 
represent small phase errors rather than relatively large fractional velocity errors. Only 15 of 
the 95 points agree to worse than 5%. Fit #2 minimizing the error with respect to individual 
kicks performs even better in this regard, with only 10 points differing by more than 5% and 
only 2 points differing by more than 10%. In either case, th e fits with the new ansatz are 



significantly better th an either the simple r rj 2 fit proposed by lLousto Zlochowerl ( 120091 ) or 
the rf fit proposed by iBaker et al.l (120081 ) . 
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Table 1. Recoil data from 



Lousto & Zlochower! (12009]) (set A), baker et all fcpQ8h (set B) 



this work (set C), iDain et al. J2008h (set D), and ICampanelli et al.l J2007bh (set E). The 
spin angles 0i and 02 are in radians, and the recoil velocities Vu are in km s -1 . The 8th 
column shows the result of a fit intended to minimize the error relative to the maximum 
recoil velocity per (q, a^, a^) triplet and the 9th column shows the relative error from that 

fit. The 10th column shows the result of a fit intended to minimize the conventional 
relative error and the 11th column shows the relative error from that fit. In both cases, the 

vast majority of points agree to well within 10%. 
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Table 1 — Continued 
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n nnn 

u .uuu 


n 78fi 

U . 1 ou 


n nnn 

u .uuu 


9 ^81 

Zj .UO X 
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Uc/U.U 


ssn 5 

uou. u 


n n9fi 

U. UZjU 
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UO 1 .o 


n nn7 

U. UU I 


R 


n sss 


n 9nn 

u.^uu 


n n99 

U . UZjZj 


n nnn 

u .uuu 


3 149 

U . XtiZj 


4Q n 

HzC/.U 


^n 5 

uu. u 


n nsi 

U. UU X 


^9 n 

UZj .u 


n nfi9 

U. UUZj 


R 


n sss 


n ?nn 

u.^uu 


n n99 

U . UZjZj 


^ 4Q8 

U .t:C/0 


9 S^fi 

Zj .UUU 


48 n 


48 1 

tiO. X 


n nns 

u. uuu 


48 Q 

tiO . C7 


n ni 8 

u. u xo 


R 


n 333 

U .tJUiJ 


n 2nn 


n n99 

U . UZj Zj 


4 71 2 


1 F571 

X . U 1 X 


1 7 n 

X l .u 


1 7 6 

X 1 . U 


n n34 

U . UUt: 


1 7 2 

X 1 .Zj 


n ni 2 

U . U X Zj 


R 


n sss 
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U.^UU 
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U . UZjZj 


n nnn 

u .UUU 


n nnn 

u .UUU 


1 1 4 n 
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n nns 

u. UUU 


1 1 ^ 9 
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R 
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-S7 n 
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u. uuu 
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UU. 1 
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u. uuo 


R 
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u .uuu 


n 9nn 

u.^uu 


n n^n 

u .uuu 
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U .t:C/0 


9 S^ifi 

Zj .UUU 


1 1 1 n 

X X X . u 


1119 

X X X . Zj 


n nn9 

U. UUZj 


114 1 

X It. X 


n n98 

U. UZjO 


R 


n ^nn 

u . uuu 


n 9nn 

U .ZjUU 


n n^n 

u .uuu 


4 71 9 


1 ^71 

X .U 1 X 


1 qs n 

X C/U.U 


1Q5 5 
X C/U. u 


n m s 

U. U X u 


1 Q8 n 

X C/O .u 


n n9fi 

U. UZjU 


R 


n ^nn 

u .uuu 


n 9nn 


n n^n 

u .uuu 


Fi 498 

U .rtU O 


1 571 

X . U 1 X 


75 n 

1 u.u 


73 2 

1 U . Zj 


n n23 

U . UZjU 


76 4 


9 91 9 


R 


n ^nn 

u .UUU 


n 9nn 

u.^uu 


n n^n 

u .UUU 


n nnn 

u .UUU 


1 571 

X .U 1 X 


-55 n 

uu .u 


-5fi 7 

UU. 1 


n nsi 

U. UU X 


n 

uu.u 


n nnn 

u. UUU 


R 


u .uuu 


n 9nn 

u^uu 


n nsn 


i n47 

X .U1 I 


4 18Q 


-S81 n 

UO X .u 


-S89 Q 

UOZj. C/ 


n nn^i 

u. uuu 


-S84 ^ 


n nnn 

U. UUC/ 


R 


u .uuu 


n 9nn 

u.^uu 


n nsn 


n nnn 

u .uuu 


S 149 

U . XtiZj 


-i s^ n 

xuu .u 


-1 S9 1 

XUZj . X 


n n99 

U . UZjZj 


-1 SS 4 

XUU.t: 


n ni 9 

U. U XZj 


R 


u .uuu 


n 9nn 

u.^uu 


n nsn 


^ 4Q8 

U .t:C/0 


9 S5fi 

Zj .UUU 


1 R8 n 

xuo.u 


1 fi5 s 

xuu. u 


n m fi 

U. U X u 


1 fi^i 9 

XUU .Zj 


n m 7 

U.Ul 1 


B 


n.666 


9.299 


9.989 


4.712 


1.571 


364.n 


365.9 


n.nns 


367.n 


n.nns 


B 


n.769 


9.299 


9.118 


9.999 


3.142 


-386.n 


-387.2 


n.nn3 


-388.1 


n.nns 


B 


n.769 


9.299 


9.118 


5.498 


2.356 


-525.n 


-521.7 


n.nn6 


-521.2 


n.nn7 


B 


n.769 


9.299 


9.118 


4.712 


1.571 


-348.n 


-359.6 


n.nns 


-349.9 


n.nn3 


B 


9.999 


9.299 


9.165 


9.999 


3.142 


-542.n 


-542.8 


n.nm 


-542.4 


n.nm 


B 


n.999 


9.299 


9.165 


5.498 


2.356 


-657.n 


-656.3 


n.nm 


-655.2 


n.nn3 


B 


n.999 


9.299 


9.165 


4.712 


1.571 


-384.n 


-385.3 


n.nn3 


-384.1 


n.nnn 


C 


9.674 


9.367 


9.236 


2.513 


5.655 


-973.7 


-944.4 


n.n3n 


-941.9 


n.n34 


C 


n.674 


9.367 


9.236 


2.994 


4.189 


-442.9 


-435.8 


9.916 


-431.8 


n.n25 


c 


n.674 


9.367 


9.236 


2.994 


9.999 


-988.5 


-934.6 


9.955 


-934.4 


n.nss 
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Table 1 — Continued 



set 


Q 


1 




01 


02 


num. V\\ 


fit Vf| #1 


1 A t r /t r 1 

|AV/V| 


fit Vf| #2 


1 A T r /t t 1 

|Avy\/| 


V 


n fi74 


n 3fi7 





236 


1 957 

X .ZjO 1 


4 3Q8 


-3fi1 7 

OU X . ( 


-344 3 

Ot:t:.0 


n n48 


-33Q 4 


n nfi9 

U. UUZj 


n 

V 


fi74 


n 3fi7 





236 


fi nnn 


4 1 8Q 

L ±. XOC/ 


374 9 


331 fi 

OO X . u 


n 1 14 


33fi 5 
oou.o 


n mi 

U. XU X 


c 


D fi74 


n 3fi7 





236 


n nnn 


3 149 

O. XtiZj 


74Q n 

I H:C/ .U 


731 fi 


n n93 


731 9 

I O X . Zj 


n n94 

U. UZjH: 


n 


fl fi74 







236 


n nnn 


9 QQ4 

Zj. Ujt: 


fi79 4 

U 1 Zj .4: 


fi89 9 

UOZj . Zj 


n m 5 

U .U XtJ 


fi77 3 

U I I . o 


n nn7 

U. UU I 


C 


fi74 


n 3fi7 





236 


5 n97 

O .UZj I 


1 885 

X . OOO 


89fi 1 

OZjU . X 


7Qfi 5 

1 C/U. O 


n n3fi 


7Q1 4 


n n49 

U. UtiZj 


n 


fi74 







236 


4 1 HQ 

L ±. XOC/ 
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fin3 n 
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Oj 1 .J 
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U. UOO 


n 


fl fi74 
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4 1 HQ 
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-945 fi 
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ZjtiU . t: 


n nn3 
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n nnn 

u. uuu 


n 

V 


n fi74 


fl 3fi7 





236 


3 77Q 
o. / / u 


n fi98 

U. UZjO 


-937 7 

ZjO I . 1 


-93Q 4 


n nn7 


-942 1 

ZjtiZj . X 


n ni q 

U. U X ZJ 


c 

V 


U .O X O 


fl 3fi7 





177 


9 51 3 

Zj .O lO 


5 fi55 

O. UOO 


-584 3 

OOti.O 


-578 n 

O 1 o.u 


n ni 1 


-587 n 

OO I .u 


n nn5 

u. uuo 


f! 


n 51 5 

U .O X O 


fl 3fi7 





177 


9 nQ4 

Zj .\JCJ l ± 


4 18Q 

L ±. XOC/ 


-i 93 n 

- XZjO.U 


-1 n3 3 


n 1 fin 


-118 1 

X XO . X 


n n4n 

U. Ut:U 


n 

V 


n 51 5 

u . o x o 


fl 3fi7 





177 


9 QQ4 

Zj .Ifjl 


n nnn 


-fi84 8 

UOt.O 


-fi83 1 

UOO. X 


n nn9 


-fi85 n 

UOO .u 


n nnn 

u. UUU 


c 


51 5 

U . O X O 


367 

W . OU 1 





177 


1 957 


4 3Q8 

t:. OJU 


-1 fi Q 

X U . U 


1 n 

X . u 


1 n5Q 


-1fi 5 

xu . o 


n n24 

U . UZjt: 


c 

V 


n 51 5 


fl 3fi7 





177 


n nnn 


4 1 8Q 

t:. XOC/ 


483 Q 


4fi9 fi 

t:UZ. U 


n n44 


448 5 


n n73 

U.Ul o 


n 

V 


n 51 5 

U .0 x o 


fl 3fi7 





177 


n nnn 


3 149 

O. llz 


57Q 8 

O I ZJ .O 


578 fi 

O 1 O.U 


n nn9 


57fi 8 

O I u. o 


n nn5 

u. uuo 


n 


n 51 5 

u .0 x o 


fl 3fi7 





177 


n nnn 


9 QQ4 

Zj. UC/t: 


345 9 

OIO .Zj 


34fi 7 


n nn4 


357 5 

OO I . o 


n n3fi 

u. uou 


n 


n 51 5 

u . o x o 


fl 3fi7 





177 


5 n97 

O .Uz 1 


1 885 

X . OOO 


3fi9 4 


35fi fi 


n mfi 


373 n 

Ol o.u 


n n3n 

u. uou 


v 


51 5 

U . O X O 


367 

W . OU 1 





177 


4 1 8Q 

t: . XO <-/ 


2 QQ4 

Zj . WJt 


2Q2 1 

Zj WZj . _1_ 


22n 5 

ZjZj w . o 


n nQi 


23fi fi 

ZjOU . u 


n 1 71 

U. X 1 X 


n 

V 


n 51 5 

U . O X O 


fl 3fi7 





177 


4 1 HQ 

L ±. XOC/ 


n nnn 


-931 1 

ZjO X . X 


-943 4 


n n53 


-93Q 4 

ZjOC/ .t: 


n n3fi 

u. UOU 


c 

V 


n 51 5 

U .O X o 


fl 3fi7 





177 


3 77n 

O.I (U 


n fi98 

U. UZjO 


-354 8 

OOti.O 


-358 9 

OOO. Zj 


n nm 


-34fi 3 

Ot:U.O 


n n94 

U. UZjH: 


n 

1 _9 


1 nnn 

X .www 


f) Q3Q 





930 


1 571 

X . O 1 X 


4 71 2 


9379 n 

ZjO 1 z.U 


941 3 fi 

Zjt: X O . U 


n ni8 

\J .u xo 


9493 7 

ZjtiZjO. 1 


n n99 

U . UZjZj 


n 

1 J 


1 nnn 

X .www 


fl Q3Q 





930 


1 833 

X . ooo 


4 Q74 


9887 n 

ZjOO I .u 


9Q79 9 


n n3n 

u.uou 


9Q75 8 

ZjC/ I o. o 


n n3i 

U. UO X 


D 


1.000 


0.930 





930 


2.269 


5.411 


3254.0 


3440.6 


0.057 


3432.8 


0.055 


D 


1.000 


0.930 





930 


3.176 


0.035 


2226.0 


2390.4 


0.074 


2366.0 


0.063 


D 


1.000 


0.930 





930 


4.817 


1.676 


-2563.0 


-2659.2 


0.038 


-2666.8 


0.041 


D 


1.000 


0.930 





930 


4.974 


1.833 


-2873.0 


-2972.2 


0.035 


-2975.8 


0.036 


D 


1.000 


0.930 





930 


5.149 


2.007 


-3193.0 


-3233.9 


0.013 


-3232.9 


0.013 


D 


1.000 


0.930 





930 


5.934 


2.793 


-2910.0 


-3152.3 


0.083 


-3133.1 


0.077 


E 


1.000 


0.515 





515 


1.571 


4.712 


1833.0 


1881.1 


0.026 


1875.6 


0.023 


E 


1.000 


0.515 





515 


0.785 


3.927 


1093.0 


1079.6 


0.012 


1076.4 


0.015 


E 


1.000 


0.515 





515 


3.142 


0.000 


352.0 


354.4 


0.007 


353.3 


0.004 


E 


1.000 


0.515 





515 


4.712 


1.571 


-1834.0 


-1881.1 


0.026 


-1875.6 


0.023 
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5. Ejection probabilities and discussion 



One of the most important outputs of kick calculations and fits is the probability dis- 
tribution of kicks given assumptions about the mass ratio, spin magnitudes, and spin direc- 
tions. Tlus_o^tnbjition is critical to studies of hierarchical merging in the early universe 
(e.g., IVolonteril ( 120071 )) as well as to the gas within galaxies ( IDevecchi et al.l 120091 ) and an 



evalu ation of the prospects for grow th of intermediate- mass black holes in globular clus- 
ters (IHolley-Bockelmann et al.ll2008l ). In Table 2 we show the resu l ts of o ur work (from fit 
#1), compared with the proposed fit formula of ICampanelli et al.l (12007a[). It is clear that 
our work gives distributions very close to those of Campanelli et al. ( j2007a ). with perhaps 
slightly smaller kicks because of the rf term we include. 

In summary, we have demonstrated that a modified formula fits all available out-of-plane 
kicks extremely well. The wide range of mass ratios, spin magnitudes, and angles explores all 
the major aspects of parameter space for the out-of-plane kicks, and thus we do not expect 
new results to deviate significantly from our formula. The excellence of these fits suggests 
that the kick distribution is known to an accuracy that is sufficient for any astrophysical 
purpose. 



New simulations used for this work were performed on Jaguar at Oak Ridge National 
Laboratory. MCM acknowledges partial support from the National Science Foundation un- 
der grant AST 06-07428 and NASA ATP grant NNX08AH29G. The work at Goddard was 
supported in part by NASA grant 06-BEFS06-19. We also wish to think S. McWilliams and 
A. Buonanno for helpful discussions. 



Table 1 — Continued 



set 


q 


a>i 


a 2 


01 


02 


num. V\\ 


fit V\\ #1 


|AV/V| 


fit V\\ #2 


\*v/v\ 


E 


1.000 


0.515 


0.515 


3.304 


0.162 


47.0 


46.3 


0.014 


46.2 


0.017 


E 


1.000 


0.515 


0.515 


0.000 


3.142 


-351.0 


-354.4 


0.010 


-353.3 


0.007 



Table 2. Fraction of k ick speeds above a given threshold, compared with the results of 
Campanelli et al.l ( ]2007ai ) (CLZM). In all cases we assume an isotropic distribution of spin 

orientations. 



Mass ratio and spin Speed threshold CLZM This work 



1/10 < 


1 < 1 


, Oil = a? = 


= 0.9 


v > 500km 


s" 


-1 


0.364±0.0048 


0.342526±0. 00019 










v > 1000km 


s" 


-1 


0.127±0.0034 


0.120974±0. 00011 


1/4 <q 


< 1, 


Oil = «2 = 


0.9 


v > 500km 


s" 


-1 


0.699 ± 0.0045 


0.697818±0.00026 










v > 1000km 


s" 


-1 


0.364±0.0046 


0.353393±0.00019 


1/4 <q 


< 1, 


< qi, a 2 


< 1 


v > 500km 


s" 


-1 


0.428±0.0045 


0.415915±0.00020 










v > 1000km 


s" 


-1 


0.142±0.0034 


0.134615±0.00012 
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